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ABSTRACT 

We present a detailed analysis of three extremely strong intervening damped Lyman-a systems (ESDLAs, with log)V(Hi) > 21.7) 
observed towards quasars with the Ultraviolet and Visual Echelle Spectrograph on the Very Large Telescope. We measure overall 
metallicities of [Zn/H] — 1.2, -1.3 and-0.7 at respectively z abs =2.34 towards SDSS J214043.02-032139.2 (logV(Hi) = 22.4 ± 0.1), 
Zabs=3.35 towards SDSSJ145646.48+160939.3 (logV(Hi) = 21.7 ± 0.1) and z abs =2.25 towards SDSSJ015445.22+193515.8 
(logV(Hi) = 21.75 ± 0.15). Iron depletion of about a factor 15 compared to volatile elements is seen in the DLA towards 
J2140-0321, while the other two show deletion typical of known DLAs. We detect Hi towards J2140-0321 (log7V(H 2 ) = 
20.13 ± 0.07) and J1456+1609 (log (V(H 2 ) = 17.10 ± 0.09) and argue for a tentative detection towards J0154+1935. 

Absorption from the excited fine-structure levels of Oi, Ci and Sin are detected in the system towards J2140-0321, that has the 
largest H i column density detected so far in an intervening DLA. This is the first detection of O i fine-structure lines in a QSO-DLA, 
that also provides us a rare possibility to study the chemical abundances of less abundant atoms like Co and Ge. Simple single phase 
photo-ionisation models fail to reproduce all the observed quantities. Instead, we suggest that the cloud has a stratified structure: H 2 
and C i likely originate from both a dense (log hh ~ 2.5 — 3) cold (80 K) and warm (250 K) phase containing a fraction of the total H i 
while a warmer ( T > 1000 K) phase probably contributes significantly to the high excitation of O i fine-structure levels. The observed 
Ci/H 2 column density ratio is surprisingly low compared to model predictions and we do not detect CO molecules: this suggests 
a possible underabundance of C by 0.7 dex compared to other alpha elements. The absorber could be a photo-dissociation region 
close to a bright star (or a star cluster) where higher temperature occurs in the illuminated region. Direct detection of on-going star 
formation through e.g. NIR emission lines in the surrounding of the gas would enable a detailed physical modelling of the system. 

Key words. Quasars: absorption-lines 


1. Introduction 


Diffuse gaseous clouds in the Universe are primarily described 
by their neutral hydrogen column density, which is a quantity 
directly and accurately measurable in absorption against back¬ 
ground sources. The study of Ly a absorption lines in quasar 
spectra currently allows us to detect H i over ten orders of mag¬ 
nitude in column density, from Al(H i) ~ 10 12 to 10 22 cm” 2 , and 
over a wide range of redshifts. 

Several observations at low redshift show that clouds_with 
logAl ( Hi) > 15 tend to cluster around galaxies (e. g. iTripp et aki 
1 19981: iBowen et al.1 l2002t iProchaska et al.l 1201 ll) . In particu¬ 
lar, Lyman limit systems (LLS), with logA 7 (Hi) ~ 17 - 18, 
likely trace cool circumgalactic environments, with a possi¬ 
ble bimodality in metallicity that would represent outflows of 
metal-rich gas and accretion of metal-poor gas by galaxies 


* Based on observations collected with the Ultraviolet and Visual 
Echelle Spectrograph on the Very Large Telescope at the European 
Organisation for Astronomical Research in the Southern Hemisphere, 
Chile, under Programme ID 091.A-0370(A). 


dLehner et al.ll2013h . An inflexion in the Hi column density dis¬ 
tribution function is se en at log ATH i) ~ 19, due to the onset of 
self-shielding dPetitiean et al.lll992l) . When the column density 
reaches log ATH i) ~ 20, the gas is mostly neutral and produces 
strong Ly a absorption with prominent Lorentzian wings. The 
corresponding absorbers are called Da mped Lyma n-a systems 
(DLAs) when log N(Hi) > 20.3 (IWolfe et al][T986l) . These sys¬ 
tems are thought to arise when the line of sight passes at small 
impact parameter to a foreground galaxy. Early evidence for this 
came mostly from the similarity with Hi column densities in 
nearby galactic discs and the presence of metals at different lev¬ 
els of chemical enrichment: the smoking gun of star formation 
activity. Interestingly, this A(H i)-threshold also corresponds to 
the value above which Gal actic gas c louds have high fraction 
of cold neutral medium dKanekar et al.ll20TTl) . a prerequisite for 
star formation. Naively, one can assume that the higher the Hi 
column density, the closer the association between the absorbing 
gas and a galaxy. 

High redshift DLAs have been intensively studied (see 
IWolfe et al.l 120051 for a review), with the main scope to put 
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constraints on the formation and evolution of galaxies in the 
early Universe as well as understanding the links between star- 
formation activity and the distribution, physical and chemi¬ 
cal state of the gas. Systematic surveys for DLAs yielded 
a measurement the cosmological density of neutra l gas (e.g. 
iProchaska & Wolfel l2009t INoterdaeme et alJ l2009bl) . which is 
a key ingredient to constrain hydrodynamical simula tions of 
galaxy formation dPontzen et alfcOOHlAltav et al.l2013l) and the 
understanding of gas consumption by star formation as well as 
feedback effects. Because metals are released in the ISM by the 
stars with differe nt time-scales, studies of both the over all metal 
abundances (e.g. iRao et all 120051: iRafelski etakll2012l) and the 
abundance ratios (e.g. Petitiean et alJ 2008 : iPettini et al.l 12008b 
IZafar et akll2014l) constrains the star-formation history. It is also 
possible to derive the physical state of the gas from the exci¬ 
tation of atomic and molecular species. In particular, the sur¬ 
rounding UV field can be estimated, thereby constraining the 
instantaneous star-fo r mation rates (e.g. INoterdaeme et alj|2007t 

004~ 


ISrianand et al.l 12005 : I Wolfe et aD 120041) . Overall, DLAs are a 
very important probe of various physical proc esses that are at 
play in the formation of galaxies (e.g. lBird et al.ll2014l) . 

A steepening of the V(H i)-frequency distribution at the very 
high column density end has long been discussed. It is possi¬ 
ble that at least a fraction of high column density systems (in 
particular those with high-metall icity) a re missed in magnitude- 
limited QSO surveys (e.g. IBoisse et al.l l l998 ). However, radio- 
selected DLA samples appe ared not to differ significantly from 
optica lly-selected ones (e.g. lEllison et al.l[2005t Ilorgenson et alJ 
120061) . at least in terms o f H i-content. Another explanation, pro¬ 
posed bv ISchavel d200 ll) . is that the high column density gas col¬ 
lapses into molecular hydrogen eventually leading to star for¬ 
mation. Somehow paradoxically, this important mechanism at 
the heart of galaxy formation is observationally supported by 
the absence of high column densities, molecular gas and in-situ 
star-formation in DLA samples (see iKrumholz et a1.1(20091) . It is 
therefore important to directly search for these three ingredients 
and study the relation between them. 

Molecular gas has a very small covering fraction, making 
it unlikely to be found alo ng a given random line of sight 
(iZwaan & ProchasLal 12006 ). Nevertheless, efficient searches 
based on the presence of strong neutral carbon lines - which 
is an indicator of cold gas - have allowed us to detect for the 
first time translucent clouds (the regime in-between dif fuse and 

Snow&McCafll 2006|) in high 


dense molecular clouds, see e. g. ISn ov 
redshift DLAs (ISrianand et al.l l2008: 


Noterdaeme et al.ll2009aL 


1 20101.1201 ll) . While the A(H ^-distribution of Ci-selected ab¬ 
sorbers is found to be flatter than that of regular DLAs (i.e. Hi- 
selected), the former population is still dominated by low-/V(H i) 
and the current sample (66 C i systems, iLedoux et al. l 120141) 
does not contain very large Hi column densities. The strong 
reddening by these systems also reinstates the possibility of a 
dust bias against the cold phase. Finally, C i systems frequently 
present the 2175 A absorption feature, although with a reduced 
strength compared to Milky Way lines of sight, possibly because 
of strong UV field originating from in-situ star-formation. 

The detection of very high column density DLAs is pos¬ 
sible through searching large datasets of q uasar spectra. Us ¬ 
ing the Sloan Digital Sk y Survey (SDSS, lYork et all l2000l) . 
INoterdaeme et al.l (I2009bl) showed that the slope of the A(Hi) 
distribution at the high end is actually shallower than previ¬ 
ously thought, with the first detection of a QSO-DLA reach- 
ing lo g A(H i) = 22 (in which H 2 is detected, iGuimaraes et al.l 
120121). The tremendous inc rease of available quasar spec¬ 
tra rtParis et al.l I20121 2014 ) in the SDSS-III Baryon Oscil¬ 


lation Spectroscopic Survey (BOSS, iDawson et al.l 120131 ) al- 
lowed to extend the stud y to even higher column densities 
(INoterdaeme et al.l 1201 2bl) . with a slope similar to what is 
seen from_opacity-corrected 21-cm maps in the local Universe 
dBraunl 120121) . The analysis of stacked spectra showed that 
extremely strong DLAs (ESDLAs, with log /V(H 1 ) > 21.7) 
are likely to arise from galaxies at small impact parameters 
(b < a few kpc) with average star forma tion rates of the or¬ 
der of 2 M 0 yr 1 ( INoterdaeme et al.ll2014l) . X-shooter observa¬ 
tions of a system with logiV(Hi) = 22.10 towards the quasar 
SDSS J113520.39-001053.5 resulted in the detection of a star- 
forming galaxy at b < lkpc fro m the quasar line o f sight 
(INoterdaeme et al.ll2012ai see also iKulkarni et akll2012t) . mak¬ 
ing it possible the first mass d eterm ination of a DLA galaxy by 
weak lensing dGrillo & Fvnboll2014l) . 

Here, we search for molecular hydrogen in extremely strong 
DLAs, using high-resolution spectroscopic observations. We 
present the observations and data reduction in Sect. [2] and the 
abundance measurements in Sect. [3] We then investigate the 
chemical and physical conditions in the highest /V(H 1 ) QSO- 
DLA detected so far in Sect. [4] We summarise and conclude in 
Sect. [5] 


2. Observations and data reduction 

The three quasars were observed in service mode with the Ul¬ 
traviolet and Visual Echelle Spectrograph (UVES. Ibekker et~akl 
l2000l) mounted on the Very Large Telescope UT2 of the Eu¬ 
ropean Southern Observatory under Programme ID 091.A- 
0370(A). The total 15 600 s of science exposure of each ob¬ 
ject was divided into four exposures of 3900 s each. All ob¬ 
servations of SDSS J214043.02-032139.2 (z e m = 2.48, here¬ 
after J2140-0321) and SDSS J015445.22+193515.8 (z em = 
2.53, hereafter J0154+1935) were performed using the stan¬ 
dard beam splitter in 390+580 setting that covers roughly 
from 330 to 450 nm on the blue CCD and from 465 to 
578 nm and 583 to 680 nm on the two red CCDs. Obser¬ 
vations of SDSS J145646.48+160939.3 (z em = 3.68, hereafter 
J1456+1609) included three exposures using the 390+760 set¬ 
ting that roughly covers wavelength ranges of 567 to 754 nm 
and 763 to 946 nm for the two red CCDs along with an expo¬ 
sure in 390+580 setting. The slit width of 1 arcsec and a CCD 
readout with 2x2 binning used for all the observations resulted 
in a spectral resolution power R « 48 000 dispersed on pixels of 
~ 2.0 km/s. The log of the observations is provided in Table [T| 
We reduced the raw frames to extract the ID spectra us¬ 
ing the UVES Common Pipeline Library (CPL) data reduction 
pipeline release 6.30 by implementing the optimal extraction. 
Polynomials of 4 th order were used to find the wavelength so¬ 
lution and the rms calibration error was found to be always 
less than 100 ms -1 . We shifted the individual science expo¬ 
sures to the heliocentric-vacuum frame correcting for the obser¬ 
vatory’s motion towards the line of sight at the j:xf>osure_mid 
point and using the air-to-vacuum relation from ICiddor (11 9961) . 
We then rebinned individual spectra to a common wavelength 
array, scaled and combined them altogether using weights es¬ 
timated from the errors in each pixel while rejecting bad pix¬ 
els and cosmic-ray impacts at the same time. In the follow¬ 
ing, abundances are expressed relative to solar, i.e. [X/H] = 
log N(X)/N(H) - log(X/ //) 0 , w here we adopted the solar abun¬ 
dances of lAsplund et al.l ( 2009 ). 


1 http://www.eso.org/sci/facilities/paranal/instruments/uves/d 
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Table 1 . Log of VLT/UVES observations 


QSO 

Observing dates 

Setting 

Exposure time 

Seeing (") 

Airmass 

SDSS J214043.02-032139.2 

15-06-2013/02-07-2013 

390+580 

4 x 3900 s 

0.7-0.8 

1.1 

SDSS J145646.48+160939.3 

06-05-2013/12-05-2013 

07-06-2013 

390+760 

390+580 

3 x 3900 s 

1 x 3900 s 

0.9 

1.0 

1.3-1.4 
1.3 

SDSS J015445.22+193515.8 

01 -09-2013/04-09-2013/05-09-2013 

390+580 

4x3900s 

0.8-1.0 

1.4 




-10000 -5000 0 5000 10000 

Relative velocity (km s -1 ) 




Relative velocity (km s ’) Relative velocity (km s ) 


Fig. 1 . Voigt-profile fitting to the damped Lya- (top) and Ly/? lines at 
z a b s = 2.34 towards J2140-0321. The normalised spectrum has been 
boxcar-smoothed by 10 pixels for presentation only. We note that here, 
and in all figures representing Voigt-profile fits to the data, the y-axis 
correspond to the normalised flux (unitless). 


3. Abundance measurements 

3.1. Neutral hydrogen 


The DLA at z a b s = 2.34 towards J 2140-0321 presents a very 
wide dark core and extended wings that decre ase the quasar 
continuum over several hundred angstroms. iNoterdaeme et al. 
(120141) used a principal component analysis (see iParis et al 


12011 ) on the low-resolution, flux-calibrated BOSS spectrum to 
get a first estimate of the unabsorbed quasar continuum, includ¬ 
ing the Ly a emission line. They derived log/V(H t) = 22.35 
by fitting the damped Ly a line. Fitting both the damped Ly a 
and Ly /3 absorption lines in the UVES spectrum together with 
the quasar continuum, we derive log /V(Hi) = 22.40 + 0.10, see 
Fig.ffl This represents the highest column density ever found in 
a high-- intervening cloud. 

Similarly, we derive the Hi column densities in the DLAs 
at Zabs = 3.35 towards J 1456+1609 and z a b s = 2.155 towards 
J 0154+1935. For J 1456+1609, we use the Lya, Ly// and Lyy 
to constrain the fit while only Lya is covered with sufficient 
signal-to-noise ratio for J 0154+1935. We obtain logA(Hi) = 
21.70 + 0.10 and 21.75 + 0.15, respectively. These values are 
also in good agreement with the respective values of log /V(H i) = 
21.85 and 21.77 derived from the BOSS data. 


3.2. Metallicity and dust depletion 

We analyse the three spectra using standard ^-minimisation of 
multiple Voigt-profiles performed with vpfit v9./B The velocity 
profiles of the metals and the corresponding best Voigt-profile 
fitting are shown in Fig. [4] For the three systems, the metal pro¬ 
files are rather compact and the decomposition into individual 
components is not unique. Therefore, we quote the total col- 


2 Carswell http: //www. ast. cam. ac. uk/~rfc/vpfit. html 


Fig. 2. Same as Fig. ^ for the DLA towards J1456+1609. Boxcar- 
smoothed by 5 pixels for presentation only. 



Relative velocity (km s 1 ) 


Fig. 3. Same as Fig. ^ for the DLA towards J0154+1935. Boxcar- 
smoothed by 10 pixels for presentation only. 


umn densities in Table [2] i.e. summed over all velocity compo¬ 
nents, which should in principle be more reliable and depend less 
on the exact velocity component decomposition. We also inde¬ 
pendently checked our results using the apparent optical dept h 
method (AOD, see ISavage & Sembachlll99Tl iFox et al.ll2005h . 
The weighted-mean values obtained by applying this technique 
to several unblended transitions are in very good agreement with 
those based on Voigt-profile fitting. 

Since Znn is a volatile species, it is expected to be unde¬ 
pleted onto dust grains and to be one of the best indicator of 
the true metallicity. We do not apply any ionisation correction, 
which should actually be very small given the very large H i col¬ 
umn densities in the clouds (see e.g. lPeroux et al.ll2007l) . 

The metallicity measurement is particularly difficult for 
the system towards J2140-0321: a two-component model with 
log V(Znn) ~ 13.4 provides a good fit to the bulk of the profile. 
However, Zn H/12026 is very strong and unfortunately the only 
zinc line covered by our spectrum. This makes the column den¬ 
sity highly uncertain, despite a satisfactory fit. A more careful 
look at the metal profile, in particular Ni n lines, reveals a central 
narrow component. We find a fitting solution that converges with 
a very narrow component (b ~ 1.2 kms' 1 ) that contains a sig¬ 
nificant amount of Znn: we then get a total log V(Znn) ~ 13.7, 
i.e. a factor of two higher than above. 

The uncertainty on the Znn column density becomes also 
much larger and we get [Zn/H]-1.24+0.93. Unfortunately the 
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AOD method cannot be used to derive N (Zn n) from Zn H/12026 
because the flux in this line reaches zero values. A higher res¬ 
olution spectrum covering also Zn n/12062 is therefore required 
to get a more precise abundance of zinc. Our spectrum covers 
the position of Pn/U1301,1532 which are fortunately not satu¬ 
rated, and from which we derive [P/H]-1.05. Bec ause p hos- 

phorus is likely a non-refractory element (see e.g. lMolaro et alJ 
1200 ll . and references therein), we can use it as an indicator of the 
metallicity for this system. For the DLAs towards J1456+1609 
and J0154+1935, we find overall metallicities of respectively 
[Zn/H] = —1.3 and -0.7 based on unsaturated Zn n/12026 and 
Zn n/12062 absorption lines. 

Tin, Cru, Nin and Fen are refractory elements that de¬ 
plete easily onto dust grains (e.g. iPettini et al .11 19971) . We ob¬ 
serve iron depleted by a factor about 15, 3 and 3 compared to 
non-depleted elements (Zn, P) in the systems towards respec¬ 
tively, J2140-0321, J1456+1609 and 10154+1935. The deple¬ 
tion patterns (Fig.0 for the systems towards J1456+1609 and 
J0154+1935 are si milar to that derived fro m stacking BOSS 
spectra of ESDLAs (iNoterdaeme et al.ll2014l). This is typical o f 
what is seen in the Small Magellanic Cloud (IWeltv et al.ll 19971) . 
which also presents high H i column densities and relatively low 
metallicities. In turn, the depletion pattern in the DLA towards 
J2140-0321 is more similar to that seen in the warm ISM in the 
Milky-Way disc. 


0.0 


-0.5 


- 1.0 


-1.5 


Fig. 5. Depletion patterns of different elements. For J2140-0321 we 
use the phosphorus abundance (relative to solar) as a proxy for that 
of zinc. The grey bars correspond to the val ue ob tai ned fro m averag - 
ing low-resolution BOSS spectra of ESDLAs ( INoterdaeme et alj20141 ). 
The dotted (resp. dashed) line corresponds to the Milky-Way Halo (resp. 
Warm disc) depletion pattern. 
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3.3. Molecular hydrogen 

Molecular hydrogen absorption lines are unambiguously de¬ 
tected in the systems towards J2140-0321 and J1456+1609 
from the Lyman (B 1 ^ <— A 1 !/) and/or Werner (C l II„ <— X '£+) 
electronic bands. In the former case, very strong Lyman-band 
absorptions are seen up to rotational level J - 5 (Fig. [6|, with 
damping wings for the first two rotational levels in all transitions 
and for several J - 2 lines. This allows for accurate column 
density measurements in these low rotational levels. In turn, all 
/ = 3 lines are on the flat part of the curve of growth, implying 
a higher measurement uncertainty. Some J > 4 lines are unsatu¬ 
rated, making again the column density measurement in princi¬ 
ple more reliable. We note that the H 2 absorption is well mod¬ 
elled by a single component but with slightly different Doppler 
parameters and redshifts for the different rotational levels. This 
indicates that the different rotational populations of H 2 may not 
be fully co-spatial with a possible velocity structure in the H 2 


Table 3. Molecular hydrogen in the z abs = 2.34 DLA towards 
J2140-0321. 


Rot. level 

£abs 

b 

[km s -1 ] 

log N 
[cm -2 ] 

J=0 

2.33990 

- 

19.84 + 0.09 

J=1 

2.33999 

- 

19.81 + 0.04 

J=2 

2.33997 

4.0+ 0.7 

17.96 + 0.14 

J=3 

2.33997 

4.9+ 0.7 

17.76 + 0.40 

J=4 

2.33996 

6.5 + 0.5 

15.88 + 0.26 

J=5 

2.33998 

5.7+ 1.8 

15.17 + 0.16 

J=6 

2.33998 


< 14.72 

total 



20.13 + 0.07 



log/= -1.97 ±0.17 


Table 4. Molecular hydrogen in the z abs = 3.35 DLA towards 
J1456+1609. 


h 2 

Rot. level 

£abs 

b 

[km s -1 ] 

log N 
[cm -2 ] 

J=0 

3.351829 

0.6+0.1 

16.09+0.26 

J=1 

5 ? 


16.75+0.05 

J=2 



16.29+0.10 

J=3 

” 


16.49+0.06 

J=4 

” 


15.78+0.16 

total 



17.10+0.09 
log/=-4.30+0.20 


gas that is not resolved into distinct components with the cur¬ 
rent data. The column densities in the individual rotational lev¬ 
els of H 2 are provided in Table [3] We measure a total H 2 col¬ 
umn density of log A/H 2 ) = 20.13 ± 0.07, from which we de¬ 
rive the overall molecular fraction log / = -2.14 ± 0.17, where 
/ = 2A(H 2 )/(2A(H 2 ) + A/Hi)). 

H 2 lines are much weaker towards J1456+1609, but still 
detected in a single component up to / = 4 thanks to the 
large oscillator strengths of the lines in the covered Lyman and 
Werner bands (Fig. |7] and Table QJ. We measure logA/H 2 ) = 
17.10 ± 0.09, which corresponds to log / = -4.30 ± 0.20. Since 
the errors in each rotational level are not independent, the total 
error is taken from the extremum values in each level. 

Unfortunately, the signal-to-noise ratio of our UVES spec¬ 
trum of JO 154+1935 is too low in the blue to firmly detect 
H 2 . This makes any ^-minimisation unreliable as the total % 2 
changes little irrespective of the assumed column density. Re¬ 
binning the UVES spectrum by 4 pixels, we can see however 
that the data is consistent with the presence of H 2 with column 
density in the range log V(H2)=17-19, assuming a thermal ex¬ 
citation of T — 100 K and a Doppler parameter b - 3 kms” 1 , 
typical of what is seen in other H 2 absorbers. Neutral carbon ab¬ 
sorption lines, in turn, are unambiguously detected in two com¬ 
ponents in this system, with N( Ci) = 10 14 cm -2 , see Fig.[3and 
Table [5] Given the first ionisation potential of carbon (11.2eV), 
the presence of strong C 1 ab sorptio n is_a reliable indicator of the 
presence of Hi in the cloud ( ISrianand et al.1120051) . 

4. The DLA towards J2140-0321: chemical and 
physical conditions 

In this section, we study the z a b s = 2.34 system towards 
J2140-0321 in more details, as it is particularly interesting for 
having the largest H 1 and H 2 column densities ever measured 
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J2140-0321, z a b s = 2.33995 J1456+1609, z abs = 3.35183 J0154+1935, z abs = 2.25155 



Fig. 4. Voigt-profile fit to the metal absorption lines in the three ESDLAs. The redshift that defines the origin of the velocity scale, indicated on 
the top, correspond to that of the cold gas (from Ci for J2140-0321, J0154+1935 and from H 2 for J1456+1609). The blue tick marks indicate the 
position of the individual velocity components in the fit. . 


Table 2. Overall column densities and abundances 


Species 

J2140-0321, z abs = 2.34 

log N [cm 2 ] 
J1456+1609, Za bs = 3.35 

J0154+1935, Zabs = 2.25 

Hi 

22.40 + 0.10 

21.70 + 0.10 

21.75 + 0.15 

Sin 


15.81 +0.05 (15.83) 
[Si/H] = -1.40 + 0.11 

16.42 ± 0.02 
[Si/H] = -0.84 + 0.15 

Pn 

14.76 + 0.08 (14.74) 
[P/H] = -1.05 ± 0.13 



Tin 

13.26 + 0.05 (13.27) 
[Ti/H] = -2.09 + 0.11 



Cru 


13.75 + 0.04(13.68) 
[Cr/H] = -2.22 + 0.11 

14.07 + 0.01 (14.09) 
[Cr/H] = -1.32 + 0.15 

Fen 

15.64 + 0.03 (15.60) 
[Fe/H] = -2.26 + 0.10 

15.39 + 0.06(15.44) 
[Fe/H] = -1.81 ±0.12 

16.09 + 0.02(16.13) 
[Fe/H] = -1.16 + 0.15 

Ni ii 

14.40 + 0.05 (14.36) 
[Ni/H] = -2.22 + 0.11 

14.09 + 0.03(14.11) 
[Ni/H] = -1.83+0.10 

14.62 + 0.02(14.67) 
[Ni/H] = -1.35+0.15 

Zn ii 

13.72 + 0.92 
[Zn/H] = -1.24 + 0.93 

12.94 + 0.04(12.87) 
[Zn/H] = -1.32 + 0.11 

13.66 + 0.02(13.75) 
[Zn/H] = -0.65 ±0.15 


Notes. Column densities derived from th e Apparent Opt ical Depth method are indicated in parenthesis (see text). Abundances are expressed with 
respect to photospheric solar values from lAsplund et all |2009). 


Table 5. Neutral carbon in the z a b s = 2.25 DLA towards JO 154+1935. 


Ci 

£abs 

log A 

b 

fine-structure level 


[cm -2 ] 

[kms -1 ] 

J=0 (g.s.) 

2.25152 

13.39 ± 0.06 

5.2 ± 0.7 


2.25163 

13.06 + 0.21 

6.3 ±2.0 

J=1 (*) 

2.25152 

13.47 ± 0.06 

5.2 ±0.7 


2.25163 

13.19 + 0.14 

6.3 ±2.0 

J=2 (»*) 

2.25152 

12.71 ±0.10 

5.2 ±0.7 


2.25163 

12.56 + 0.21 

6.3 ±2.0 

total 


13.95 ± 0.05 



in a high-z QSO-DLA. In addition, the presence of chlorine, as 


well as fine-structure lines of C i, O i and Si n allows us to probe 
the chemical and physical conditions prevailing in the absorb¬ 
ing gas as well as the multicomponent nature of the absorber. In 
particular, we apply standard procedures used in Galactic ISM 
studies to extract physical parameters like density, temperature 
and radiation field and compare them with those derived for the 
EE clouds in different environments. We also consider photo¬ 
ionisation models using the code cloudy to confirm our findings. 

4. 1. Heavy elements 

Different elements are released by stars through several chan¬ 
nels, depending on the stellar masses and life-times. The very 
high column densities towards J2140-0321 provide a unique op- 
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Fig. 6. Portions of the UVES spectrum of J2140-0321 covering Hi absorption lines at z abs = 2.34. The data has been boxcar-smoothed by 3-pixels 
for the sake of visual clarity only. Each wavelength region covers a Lyman band, indicated in the top left corner of each panel. The numbers in 
blue above each absorption line indicates the rotational level to which the line belongs. 


portunity to detect rare heavy elements, bringing additional clues 
to the chemical enrichment history of the cloud. 

Absorption lines from Ge n and Co n are probably detected 
in our spectrum, albeit at less than 3cr significance level and 
from single lines (FigfTOl). We estimate the corresponding col¬ 
umn densities to be logA(Gen) ~ 12.25 and logfV(Coii) ~ 
13.05 from fitting these lines with the redshifts and Doppler 
parameters fixed to those derived from the other metals. We 
then derive [Ge/H]~ -1.8 and [Co/H]~ -2.3. Germanium is 
the heaviest element detected in our s pectrum, and has very 
rarely been d etect ed in DLAs to date dProchaska et al.l 120031 : 
iDessauges-Zavadskv et al .1120061) . Measuring the abundance of 
germanium is of particular interest because its nucleosynthesis 
resides at the transition between explosive synthesis of iron-peak 
elements and neutron-capture synthesis of heavier elements. In 
metal-poor stars, the_abundance of germanium seems to follow 
well that of iron (ICowan et al.l 120051) , favouring an explosive 
origin of this element at low met allicity, while the increase of 
[Ge/Fe] with metallicity (lRoederetil2012 !) indicates the onset of 
s-process production of germanium at high metallicity. Here, the 
observed ratio is super-solar (TGe/Fe l ~ 0.46), similar to what 
is seen in other DLAs (IDessauges-Zavadskv et al.ll2006[) . How¬ 
ever, iron is a refractory element that very easily depletes onto 
dust grains, making it difficult to interpret the observed ratio in 


terms of nucleosynthesis. Unfortunately, the abundance of un¬ 
depleted elements zinc and sulphur are not well constrained in 
the DLA towards J2140-0321 while Ni is also depleted. A pos¬ 
sible reference element could be phosphorus, which seems to 
be little depleted. However, its exact nucleosynthesis is still not 
fully understood dMolaro et al.l 1200 It). The observed sub-solar 
[Ge/P] ratio ([Ge/P] —0.75) could partly be due to depletion of 
germanium but could also reflect an enrichment dominated by 
massive stars. Cobalt is also rarely seen in DLAs, with the first 
detection reported by iEllison et al.l (1200 ll) . As argued by these 
authors, the similar depletion of cobalt and iron should make 
the [Co/Fe] ratio little dependent on the actual amount of dust. 
Here, we measure a solar cobalt-to-iron ratio, sim ilar to what is 
seen in the DLA tow ards Q 0948+43 ( Rao et al1l2005h . and in 
the Galactic disc (e.g. lProchaska et al.1120001) . The DLA towards 
J2140-0321 therefore appears to be an excellent place to study 
the nucleosynthesis processes at the origin of the chemical en¬ 
richment, although further observations, with higher S/N ratio, 
better wavelength coverage (in particular for Cu n that falls in a 
gap of our spectrum and for Zn n, for which only one transition is 
covered) and higher spectral resolution are required to put strong 
constraint on the abundances. 
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P. Noterdaeme et al.: VLT/UVES observations of extremely strong intervening damped Lyman -a systems 



4400 4410 4420 4430 

Observed wavelength (A) 


Fig. 7. Portion of the J1456+1609 spectrum covering H 2 absorption lines at z abs = 3.35. The number in blue above each line correspond to its 
rotational level. The Lyman and Werner bands are indicated as blue label below each J = 0 LL line. Other absorption lines from the same system 
are labelled in green (Lyy Ly 6 and Sin). 


4.2. Neutral chlorine 

Chlorine is a peculiar species whose ionisation potential (12.97 
eV) should make Cl n its dominant form in atomic gas. How¬ 
ever, in the presence of H 2 , chemical reactions make Cli the 
dominant form of chlorine dJural 119741) : Cl° results from the 
exothe rmic reaction Cl + +H7 —>HC1 + that has a very high rate 
(iFehsenfeld & Fergusonlfl974h . Several channels then allow the 
quick release of Cl 1 , in particular from collisions with H 2 or 
electrons. 


In addition, Cli is little depleted onto dust grains and 
should be an ex cellent tracer of clouds with high molecu- 


lar fractions (e. 

2 . Jura & York 19781 Sonnentrucker et al.ll2002t 

Balashev et all 

20151). Cl 1 has been detected in the translu- 


cent_cloud at Zahs = 2. 7 towards the quasar SDSS J1237+0647 
(iNoterdaeme et al.il2010t) . that has high molecular fraction fm > 
0.3 and also harbours CO molecules. Here, we detect a strong 
Cl 1/11347 line (see Fig. ITU, from which we derive logA(Cli) = 
13.61 +0.10 cm -2 when all parameters are left free. However, 
the data is noisy and that b may be overestimated. Since Cl 1 is 
expected to be related to the cold gas phase, we rather expect 
a velocity broadening more similar to that of H 2 (or C 1 ) and a 
good velocity match between these species. We therefore per¬ 


form a second fit, fixing b to 5kms _1 . This gives log N(C\ 1 ) = 
13.69 + 0.13 cm -2 . We therefore adopt the value logA(Cli) = 
13.67 ± 0.15, where the error conservatively represent the full 
range allowed by the two fits and their associated uncertainties. 
Because Cli is little depleted into dust and arises in the H 2 - 
bearing gas, its column density can provide a rough idea of the 
amount of H 1 present in that cold gas, provided its metallicity is 
known. Assuming a value equal to the overall metallicity, we es¬ 
timate that the H 2 -bearing gas contains a few times 10 21 cm~ 2 of 
the total hydrogen column density (10 22 4 cm 2 ). The molecular 
fraction in the remaining gas is likely too low to make chemical 
reactions with H 2 counterbalance the Cl 1 photoionisation rate. 

4.3. Carbon monoxide 

We do not detect CO absorption lines. In order to get a strin¬ 
gent limit, we stacked the covered COAX absorption bands 
(from (0-0) to (3-0)), weighting each band by //cr 2 , where / 
is the oscillator strength and cr c the noise in the continuum. 
This way, the signal-to-noise ratio in the stack is maximisecfl 

3 The signal here refers to the amplitude of the absorption, not the 
continuum. 
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3600 3605 3610 3615 3620 3625 3630 

Observed wavelength (A) 


Fig. 8. Portion of the UVES spectrum of JO 154+1935 covering the lo¬ 
cation of Lyman-bands of H 2 (from L(0-0) to L(5-0)). The spectrum 
has been rebinned by 4 pixel. The red synthetic profile corresponds to 
MTL) = 10 18 cnT 2 with a single excitation temperature of 100 K. The 
shaded region corresponds to ± 1 dex range in N( H 2 ) around this column 
density. 




5384 


5386 5388 = 5390 

Observed wavelength (A) 


5392 


Fig. 9. C 1 absorption lines in the z a b s = 2.25 DLA towards JO 154+1935. 


We assume an excitation temperature corresponding to radiative 
equilibrium with the cosmic micro wave bac kground, as seen in 
high-z CO absorption systems (iNoterdaeme et al.ll201 il l and de¬ 
rive N( CO) < 5 x 10 13 cm~ 2 , see Fig. HU This corresponds to 
CO/H 2 ~ 0.4 x 10 6 , more than an order of magnitude less than 
what has been measured in C O-bearing DLAs dSrianand et all 
120081 INoterdaeme et al.ll2010l) . but consistent with the ratio ob¬ 
served towards G alactic stars w ith the same FF column density 
(e.g. Fig. 4 of Burgh et al.ll2007l) . This shows that high FF col¬ 



-100 -50 0 50 100 

Relative velocity (km s -1 ) 


Fig. 10. Possible detection of Con and Gen in the DLA towards 
J2140-0321. 



-100 -50 0 50 100 

Relative velocity (km s’ 1 ) 


Fig. 11. Neutral chlorine in the DLA at z a b s = 2.34 towards 
J2140-0321. The dashed profile corresponds to a fit will all parame¬ 
ters free, while b is fixed to that of Ci for the solid profile. 



Fig. 12. Stack of covered CO AX absorption bands (from (0-0) to (3-0)) 
towards J2140. The red (resp. blue) profiles correspond to the expected 
profile having 1 time (resp. 3 times) the area sustained by noise (hori¬ 
zontal dotted lines) over the velocity range covered by the band (vertical 
lines) and assuming T tx = Tcmb- 


umn den sity is n ot a su fficient condition for the presence of CO 
(see also lBolatto et al.ll2013l) . 


4.4. Excitation of H 2 : Gas kinetic temperature and UV 
pumping 

At logATFF) > 16, the low rotational level populations are 
most likely to be in thermal equilibrium. Therefore, the cor¬ 
responding excitation temperatu re provides a good estimate of 
the g as kinetic temperature (e.g. iRov et akll2006l : iLe Petit et alJ 
120061) . In Fig. |T3] we show the population of the different ro¬ 
tational levels of FF as a function of the energy in these lev¬ 
els. The .1=0 to J=2 levels can be explained by a single exci¬ 
tation temperature of T ex ~ 75K, which is lower than what is 
generally observ ed in both high-’ and low-z FF-bearing DL As 
«7oi> ~ 150 K. ISrlanand et al.ll200l iMuzahid et akllml or 
at high Galactic latitude ((Toi) ~ 125 K, Gillmon et al.l 2006h 
but very similar to what is observed in the G alactic disc or 
in the Magellanic clouds (Tm ~ 80 K, ISavage etakl 119771 : 
iTumlinson et al.ll2002HWeltv et al.ll2012l) . The higher excitation 
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P. Noterdaeme et al.: VLT/UVES observations of extremely strong intervening damped Lyman-or systems 
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Fig. 13. Excitation diagram of H 2 towards J2140-0321. The straight 
blue (resp. green) dotted line corresponds to the excitation of the low 
(resp. high) rotational levels. The orange curve correspond to the sum 
of the two populations in Boltzmann equilibrium at the two different 
temperatures. 


temperature (T ex ~ 240 K) inferred for J>3 rotational levels 
shows that they are not thermalised if arising from the same re¬ 
gion and that an additional excitation process is at play. Gener¬ 
ally, such a situation is considered to be due to UV pumping, 
although ^-formation on d ust grains co uld contribute to pop¬ 
ulate the high energy levels (IWagenblastll 19921) . Another possi¬ 
bility could be that the high-J levels are also populated by colli- 
sional excitation, but in a warmer, presumably external layer of 
the cloud. The observed excitation diagram would then be the 
sum of a low temperature component (having a temperature of 
~75 K), contributing to most of the low-J column density and a 
warmer component (having a temperature of ~240 K), responsi¬ 
ble for the high-J column densities. If we assume such an excita¬ 
tion temperature of the warm gas then the measured column den¬ 
sities at J — 4 and / = 5 are consistent with a total log N(H 2 )~ 
19 in this component. 


4.5. Neutral carbon: volumic density of hydrogen and 
electrons in the cold phase 

Absorption lines from neutral carbon in the ground state ( J = 0, 
C i), first (J - 1, C T) and second (/ = 2, C T*) excited levels are 
detected in our spectrum (Fig. M- From fitting the absorption 
lines around 1560 and 1656 A (in the DLA’s rest-frame), that are 
located outside the Ly or-forest, we derive the total (ground and 
excited states) C i column density to be log N( C i) = 13.57 +0.03. 
Both the Ci*/Ci and CT*/Ci ratios indi cate high particle den¬ 
sity in the cloud. From Fig. 2 of ISilva & Viegasl (120021) and 
Fig. 12 of ISrianand et al.l ( 2005 ). we can already see that the 
density in the C i-bearing gas should be of the order of /; H ~ 
200crrT 3 . We note that at such density, the contribution from 
radiative excitation (UV pumping and cosmic microwave back¬ 
ground) becomes negligible compared to collisions. While Ci 
and Flo are not necessarily fully co-spatial, to the first order, we 




Observed wavelength (A) 


Fig. 14. Neutral carbon in the DLA at z abs = 2.34 towards J2140-0321. 


Table 6. Additional species at z a b s = 2.34 towards J2140-0321. 


Species 

log N 
[cm" 2 ] 

b 

[km s" 1 ] 

Ci,J=0 (g.s) 

13.03 + 0.04 

5.0 ±0.8 

Ct,J=l (*) 

13.20 + 0.04 

” 

Ci, J=2 (**) 

13.02 + 0.05 

” 

total 

13.57 + 0.03 


Cl 1 

13.67 + 0.15 

5-10" 

CO 

< 13.73 


Ot,J=2 (g.s.) 

17.9" 

11 

Oi, J=1 (*) 

13.89 + 0.11 

9.5 ± 1.0" 

Oi, J=0 (**) 

13.85 + 0.08 

9.5 ± 1.0 

Sin, J= 1/2 (g.s.) 

> 16.16 


Sin, J=3/2 (*) 

12.81+0.03 

9.5 ± 0.8 


Notes. ,a) The two values correspond to b fixed to that of C i and that 
obtained when left as free parameter. <w Value from fitting the line, con¬ 
sistent with the [Zn/H] and [P/H] ratios. See text. (c) Because OTT1304 
is noisy, we fixed b to the value obtained from fitting Oi**/I1306 alone. 


can consider this density to be that of the molecular gas. We note 
that this implies that the cloud is of pc-scale, when considering 
about 10% of Hi in this phase (Sect. I4.3l and l4~9l) . 

For the inferred metallicity, if we ignore the depletion of C 
into dust grains and a possible under-abundance of carbon com¬ 
pared to other elements, we expect N (C n) - 6 x 10 17 cm" 2 . This 
implie s V(Ci)/JV(Cn)=10" 4 2 . Using Eq. 5 of ISrianand et al.l 
( 2005 ) we derive n e /T — 1.4xl0 6 when we assume T~100 K. 
Here, T is the photo-io nisatio n rate of C°. If we take T ~ 2xlO" 10 
found in our galaxy (IWeltv et al.l2003h then we get n. - 3xl0" 4 . 
This value is much lower than the inferred range (0.7-4.7) xlO" 2 
in other high-z H i abs orbers using similar approach (see Ta¬ 
ble 4 of ISrianand et al.l 12005 ). We note however that the lower 
inferred electron density could also result from a much higher 
actual photo-ionisation rate than the Galactic mean value or the 
inferred log N(C n) is overestimated. Indeed, carbon could actu¬ 
ally be under-abundant, similarly towhat ha s be en suggested for 
the lack of C i in a low-' H 2 DLA (ISrianand et al.l2014 ). We will 
come back to this while discussing photo-ionisation models. 

4.6. Fine-structure levels of Oi 

We detect very strong 01/11302 absorption together with the ab¬ 
sorption from the two fine-structure levels of the ground state 
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Fig. 15. Neutral oxygen fine-structure lines in the DLA at z abs = 2.34 to¬ 
wards J2140-0321. The absorption lines seen at A ~ 4348 and 4356.5 A 
correspond to Pn41301 and Siii/11304, respectively. 


(Oi*41304 and O i**4 1306, see Fig.[B)0. This is, to our knowl¬ 
edge, the first detection of excited neutral oxygen in an inter¬ 
vening absorption system. However, OT and Oi** are more 
commonly detected in DLAs associated to GRB afterglow s 
dVreeswiik et al.1 l2004t IChen et all l2005t iFvnbo et~akl 1 20061) . 
where it is generally interpreted as an indicator of high volu- 
mic density. Unfortunately, the weaker 01/11355 line is not cov¬ 
ered by our spectrum and other O i lines are far in the blue and 
strongly blended in the Ly a forest. Fitting the 01/11302 absorp¬ 
tion line with a single component gives log N(Oi) = 17.9 ± 0.2. 
The Doppler parameter is left free and we find b ^ 11 l< m s 1 , 
surprisingly close to the b ~ 12 kms -1 expected from the ve¬ 
locity width of the profile (A,, ^ 28 k ms 1 . Sect. I4.81l. using th e 
relation between these quantities from lNoterdaeme et all ( 2014 ). 
The derived oxygen abundance would then be [O/H] ~ -1.2, i.e. 
consistent with that derived from zinc and phosphorus. 

We derive logA(Of) = 13.89 ± 0.11 and logA(Oi**) = 
13.85 ± 0.07 from fitting OTT1304 and Oi**41306. We note 
however that the former line is noisy and possibly affected by 
blends. We therefore fixed the Doppler parameter to the value 
from 01**41306 alone but note that this has little influence on 
the derived column density. 

Using the Hi collisional rate constant from 


Abrahamssonrtak (120071) and radiative d ecay rates from 


Galavis et al.l ( 199% . Jenkins & Tried (1201 lh have found that 
A(Oi*)/A(Oi**) can be a good indicator of local kinetic tem¬ 
perature (see their Fig 15). For T~80 K as suggested by Tq\ we 
expect the ratio to be higher than 2. However, the observed ratio 
(~ 1) is rather consistent with a gas temp erature of 1000 K. It is 
clear from Fig 6 of lSilva & Viegasl (120021) that the observed ratio 
of N(Oi*)/N(Oi) and N(Oi**)/N(Oi) cannot be reproduced 
by a single density for T - 150 K. However, this ratio can be 
reproduced for T = _1000K and n H ~ 300 cm -3 . It is also 
clear from Fig 6 of lSilva & Viegasl (120021) that the UV pumping 
may not be very efficient for the observed range in ratios of 
column densities. Therefore, the observed fine-structure level 
populations of O i favours the presence of warm neutral gas in 
addition to the dense cold gas probed by the low rotational level 
population of H 2 . Direct excitation by infra-red photons is also 
a possibility. 


4 . 7 . Excited C11 and Si 11 

Strong and saturated C if 41335 absorption is detected in this 
system, with an equivalent width of roughly 0.2 A. Because 
the Cn* column density is directly related to the gas cooling 
through C 11 158 pm emission CP 3/2 (Cn*) —> 2 Pi /2 (Cn)), it 
allows in principle to derive the surrounding UV flux by equat¬ 

4 Note that the ground state of neutral oxygen is comprised of the 
2s 2 2p 4 3 PjA ! o fine-structure levels, the ground-state of which corre¬ 
spond to / = 2 
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Fig. 16. Sin* at z abs = 2.34 towards J2140-0321. The absorption at 
+50 km s -1 in the Sin*41264 panel is due to Sin*41265. 


ing this cooling rate with the heating from photo-electric effect 
onto dust grains (IWolfe et al.li2003l) . although cosmi c-ray heat¬ 
ing may contribute significantly as well dDutta et al.ll2014l) . Un¬ 
fortunately, the strong saturation of C if 41335 prevents mean¬ 
ingful measurement or useful limit on N (C if) and the weaker 
C ii* 41037 is blended with Lya forest absorption. Nonetheless, 
the very strong absorption should still indicate significant cool¬ 
ing and hence star-formation activity. 

The excited fine-structure level of ionised silicon 
(Si n(J =3/2) or Si if) is also detected in this system. The 
first detection of this kind in an intervening system was re¬ 
ported by iKulkarni et al.l (120121 ) in the extremely strong DLA 
towards J1135—0010, which arises from a young low- impact 
param eter galaxy with SFR ~ 25M n yr -1 (iNoterdaeme et aD 
2012 a). More recently, iNeeleman et al.l d2015l) suggested that 
this line arises from highly turbulent ISM in young star¬ 
forming galaxies. Here, we unambiguously detect Sin*41264, 
Sin*41265 and Si if41533 (Fig. FFTTI) from which we derive 
logA(Siif) = 12.81 + 0.03. Assuming an intrinsic solar abun¬ 
dance ratio of silicon with respect to zinc and 0.2 dex depletion 
onto dust as seen in different galactic and DLA environments 
(Fig.0, we estimate, in the gas phase, A(Si n) to be about a few 
times 10 16 cm -2 , implying ATSi n*)/ATSi 11 ) ~ 10 4 . It is well 
known that the electron collisional cross-section is three orders 
of magnitude higher than that of hydrogen for the fine-structure 
excitation of Si n. If we consider only the collisions by electrons 
followed by radiative de-excitation we get n e ~ 0.8 cm -3 if 
Sin* absorption is associated to the H 2 gas with T-100 K 
(see lSrianand & Petitieanl1200(1) . We get n e ~ 0.06 cm -3 if we 
assume the gas temperature to be 1000 K. 

4.8. Kinematics 

The study of the velocity extent of different species can pro¬ 
vide additional view on the overall puzzle (Fig. fTTb . The top 
panel presents the measurement of the velocity width defined 
as the interval comprised between 5% and 95% of the to¬ 
tal optical depth of a low-ionisation species. We find 6v = 
28 kms -1 which is in the low range of the distribution for the 
overall populationjtf DLAs. This may be indicative of a low- 
mass system dLedouxet al.l 12006]) . The widt!0 of the Civ pro¬ 
file (FWH M~40 km s -1 ) is also among the lowest measured in 
DLAs dFox et al.ll2007 ) which suggests a quiet environment: its 
origin could be due to photo-ionisation by nearby massive stars 
rather than collisional ionisation in a hot halo or in galactic winds 

5 Because C iv is saturated, the kinematics cannot be quantified using 
the 5-95% range of the optical depth. 
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P. Noterdaeme et al.: VLT/UVES observations of extremely strong intervening damped Lyman-a systems 
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Fig. 17. Measurement of metal-absorption kinematics (using 
Fe H/U611, top panel) at z abs = 2.34 towards J2140-0321. The other 
panels represent the normalised spectrum around several species on the 
same velocity scale. 

(e.g. lOppenheimer & Davel2006l) . We may anticipate that in-situ 
star formation could contribute to the high excitation of Si n and 
Oi. These excited species also have b ~ 10 kms -1 , twice larger 
than what is seen in the cold phase {b(Ci) ~ 5 kms -1 ), hence 
possibly arising from a warmer and more turbulent gas. In con¬ 
clusion, we are likely observing cold gas in a star-forming region 
of a possibly low-mass galaxy. The presence of nearby massive 
stars would result in warming and exciting the outer layers of the 
cloud, giving rise to the observed Si iT and O T. 

4.9. Comparisons with models 

In this section, we intent to understand the abundances and ex¬ 
citation of the different species in a single picture by comparing 
the observations with models using the spectral simulation code 


Cloupy ( Ferland et al .1119981) . including the micro-physics of Hi 
(IShaw et al 120051) . 

To start with, we consider the case of absorbing gas not as¬ 
sociated to any star-forming region. In this case we assume the 
ionising radiation to be the meta-galactic UV background domi¬ 
nated by QSOs and Galaxies (see for example lHaardt & Madaul 
120121 ). We assume the metallicity and dust depletion to be identi¬ 
cal to the observed values. For simplicity we assume the gas tem¬ 
perature to be 80 K as inferred from 7 01 and did not ask the code 
to self-consistently obtain this temperature. The calculations are 
stopped when the total H 2 column density in the model reaches 
the observed value. We run the model for a range of hydrogen 
densities. These “cold” models correspond to the solid black 
lines in Fig. [18] We find that the observed abundance of chlorine 
in the cold gas can be reproduced for 1.4 < log n h < 1.8. How¬ 
ever, this density range is inconsistent with the C 1 fine-structure 
excitation and the density is also too low to explain the high exci¬ 
tation of oxygen fine-structure levels. Conversely, a higher den¬ 
sity model (log«H ~ 2.7), consistent with the Ci fine-structure 
excitation under-predicts the abundance of chlorine. This hap¬ 
pens because Hi is then formed more efficiently and reaches the 
observed column density at a much lower total gas column. Fi¬ 
nally, we note that the predicted C 1 column density in the cold 
phase is about a factor 20 above the observed values. 

Therefore, we considered another set of models where we 
also include the local radiation field. For this we considered the 
Galactic UV background given by iHabind d!968l) scaled by x- 
For simplicity, we also scale the cosmic ray density by the same 
amount. We considered^- = 1,10 and 30 (respectively dotted, 
dashed and dash-dotted curves). We also consider a set of model 
with the addition of a warmer component (T ~ 300 K) that con¬ 
tains logV(H 2 ) ~ 19, as discussed in Sect. 14.41 These models 
(“cold+warm”) are shown as orange curves in Fig. [18] 

We can see that irrespective of the radiation field and the 
inclusion or not of a warm envelope, the fine-structure excita¬ 
tion of Ci indicates a high density as derived above. The inclu¬ 
sion of a warm phase makes the different ratios in better agree¬ 
ment with each other, although C T* tends require slightly higher 
excitation than Ci*. As we already consider excitation by cos¬ 
mic microwave background, collisions and UV pumping self- 
consistently in cloudy, the additional pumping could be related 
to local infra-red radiation field not explicitly included in our 
models. However, if we allow for 1.5<x uncertainties in the mea¬ 
sured ratios then our models can consistently reproduce the ob¬ 
servations when 2.5< log ne <3.0. For this range in density the 
observed /V(C1 1 ) requires x close to or slightly below 10. In that 
case, the cold phase should contain between 10% and 50% of the 
total hydrogen column density while these values are between 
30% and 80% when adding the warmer phase. 

The addition of the warm phase improves significantly the 
agreement between the predicted and observed ratio of 0 1 fine- 
structure levels, in particular the O T/0 1 ratio. In turn. Si iT/Si 11 
becomes now over-predicted for the above x and n while the 
changes are still not good enough to explain the high O T*/0 T 
ratio. In the case of silicon, we remind that we were unable to 
measure /V(Si 11 ) directly due to line saturation, so the disagree¬ 
ment remains reasonable. We also note that a lower density in 
the warm phase (log n ~ 2) than in the cold phase would help 
making the agreement between model and observations better. 
Indeed, as noted before from Fig 6 of lSilva & Viegasl (120021) . it 
is clear that when T= 150 K (close to what we use in our mod¬ 
els) the very similar value of 10 4 found for N(Of)/N(Oi) and 
V(OT*)/V(Oi) cannot be produced by a single density. In turn, 
from the bottom panels of the same figure we can see it will 
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Fig. 18. Results of photo-ionisation models. For each hydrogen density the calculation is stopped when the observed (V(H 2 ) is reached. The blue 
shaded regions give the 1 and 2a range on the total column density (top) or column density ratios (bottom) from the observations. The solid, 
dotted, dashed and dashed-dotted lines correspond to an in-situ radiation field of respectively^ = 0,1,10, 30 times the local Galactic field. The 
cold-gas-only model (T = 80 K, logAifFF) = 20.13) is represented by black lines while the orange lines represent the model with an additional 
warm component (T = 300 K and log 1V(H 2 ) = 19). 


be easier to find a solution when the gas temperature is 1000 K. 
Therefore, it is also possible that the fine-structure excitation of 
Oi originates from an even warmer (i.e > 1000 K) phase than 
the one probed by H 2 absorption, presumably with a lower den¬ 
sity. This phase would not contribute to the observed Cl 1 column 
density but would contain the remaining H 1 . 

However, the addition of an even warmer phase without H 2 
that could explain 0 1 * and O C* will not solve the problem with 
C 1 . Also, as can be seen from the corresponding panel in Fig.[T8j 
increasing further the ionisation rate of carbon is not an option 
as it also alters the H 2 equilibrium. Reducing the recombination 
rate can be achieved by reducing the electron density that will not 
alter the H 2 equilibrium abundance. Because H 1 is self-shielded, 
the electron density in the H 2 gas will originate from the photo¬ 
ionisation of metals, grains and primary and secondary ionisa¬ 
tion by cosmic rays. Grain recombination also controls the elec¬ 
tron density. One important step towards understanding this is¬ 
sue is to measure N( C 11 ) accurately. The expected value suggests 
that C n/11334 should have detectable damping wings. Unfortu¬ 
nately, the S/N of our spectrum is still too low to firmly detect 
the damping wings and constrain N(C 11 ). Another (simpler) pos¬ 
sibility is that the gas phase abundance of carbon is much lower 


than assumed. Assuming 0.4 dex depletion of C onto dust grains 
as seen in the cold ISM, the true abundance of carbon in the 
DLA should be of the order of [C/O] —0.7 or even less for the 
gas phase column density to be consistent with the models. This 
is indeed similar to what is seen in Milky- way halo stars w ith 
metallicity similar to the present DLA (iAkerman et al.lf20Q4h . In 
addition, while [C/H] is difficult to measure in DLAs, available 
measurements tow ards very low metallicity DLAs do not rule 
out this possibility (ICooke et al.ll201 lllDutta et al .1120141) . 


5. Conclusion 

We have presented high resolution spectroscopic observations of 
three ESDLAs with log A'(H 1 ) > 21.7. We have firmly detected 
H 2 in two of them, towards J2140-0321 and 11456+1609. H 2 
is also likely present with a high column density in the third 
system, towards 10154+1935, but additional data is required to 
firmly establish this. This recommends a systematic survey for 
H 2 in EDLAs as the H 2 detection rate may be much higher even 
though the measured overall molecular fractions are not as high 
as those expected in translucent and dense molecular regions. 
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We studied in details the system towards J2140-0321, which 
has the highest H i and H 2 column densities observed till now in a 
high-/ intervening system (log A'(H 1 ) = 22.4) and presents sev¬ 
eral excited fine-structure transitions that help constraining the 
physical conditions. Cobalt and germanium are likely detected 
in our spectrum, although more data would be required to use 
these species to constraint accurately the nucleosynthesis history 
of this DLA. The abundance and excitation of different species 
shows that PC is predominantly originating from a cold dense 
phase (T ~ 80 K, log«H ~ 2.5-3) where chlorine is found in neu¬ 
tral form due to chemical reactions with PL. However, models 
indicate that the addition of warmer phases is required to better 
explain the high excitation of oxygen and silicon, consistently 
with the wider velocity profile of these species. Nevertheless, 
even considering multiple phases, all models overproduce C 1 , if 
we assume a solar abundance of carbon with respect to other el¬ 
ements. This may in turn indicate that the [C/O] ratio is actually 
much lower than solar in this DLA, similarly to what has been 
seen in Milky-Way halo stars. We conclude that the absorbing 
gas is likely made of multiple phases having temperature in the 
range 80 to 1000 K irradiated by an local radiation field due to in 
situ star formation. This is also consistent with the picture that 
extremely strong DLAs ari se at very small impact parameters 
from star-forming galaxies (INoterdaeme et al.ll2014l) . i.e. likely 
arise from their disc or inner regions. 

Finally, direct detection of star-formation activity through 
the detection of NIR emission lines should bring invaluable clues 
to complete the picture: it should allow us to better understand 
the excitation of the absorbing gas and link its properties with 
the overall properties of the DLA galaxy. 
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